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Abstract 

A representation of spatial infinity based on the properties of conformal geodesies is used 
to obtain asymptotic expansions of the gravitational field near the region where null infinity 
touches spatial infinity. These expansions show that generic time symmetric initial data with 
an analytic conformal metric at infinity will give rise to developments with a certain type 
of logarithmic singularities at the points where null infinity and spatial infinity meet. These 
logarithmic singularities produce a non-smooth null infinity. The sources of the logarithmic 
singularities are traced back down to the initial data. It is shown that if the parts of the 
initial data responsible for the non-regular behaviour of the solutions are not present, then 
the initial data is static to a certain order. On the basis of these results it is conjectured 
that the only time symmetric data sets with developments having a smooth null infinity 
are those which are static in a neighbourhood of infinity. This conjecture generalises a 
previous conjecture regarding time symmetric, conformally fiat data. The relation of these 
conjectures to Penrose's proposal for the description of the asymptotic gravitational field of 
isolated bodies is discussed. 

PACS numbers: 04.20.Ha, 04.20.Ex. 

1 Introduction 



A central issue in the general relativistic theory of isolated systems is that of the verification of 
the so-called Penrose proposal |22) . Following Friedrich [161 117j , the proposal can be formulated 
as: 

Proposal (Penrose, 1965). Far fields of isolated systems behave like asymptotically simple 
spacetimes in the sense that they can be smoothly extended to null infinity after a suitable confor- 
mal rescaling. 

That this proposal is not empty — i.e. that there are examples of radiative spacetimes — is 
known thanks to work by Chrusciel & Delay The idea behind their result was to combine 
a modification of the gluing construction of Corvino 0, which yields initial data which are 
Schwarzschildean outside a compact set, with the semiglobal hyperboloidal existence result of 
Friedrich jj^l- The resulting spacetime is then — see figure ^ — Schwarzschildean in a region, N, 
of spacetime "near null and spatial infinity" . However, at later times radiation could certainly 
be registered at null infinity. As emphasized by Friedrich, the radiation content of this spacetime 
is rather special: the Schwarzschildean nature of the initial data implies that Newman-Penrose 
constants of the spacetime are zero; on the other hand, the Newman-Penrose constants are know 
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Figure 1: Schematic representation of the asymptotically simple spacetimes constructed by 
Chrusciel & Delay. The initial data is Schwarzschildean outside a compact set. 

to be equal to the value of the rescaled Weyl tensor at timelike infinity |21[ I2U) . Hence, the 
Weyl tensor at vanishes for the Chrusciel-Delay spacetimes, thus implying a fast decaying 
gravitational field at very late times. 

In the light of the results by Chrusciel & Delay it is natural to ask how general can be the 
behaviour of initial data in a neighbourhood of spatial infinity if one is to obtain asymptotically 
simple solutions. A first insight to this point was provided by the analysis of Friedrich given in 
|15| . In order to keep the discussions in a reasonable complexity level, his analysis was restricted 
to time symmetric initial data with an analytic compactification at infinity. He found that a 
necessary condition in order to attain spacetimes that are smoothly extendible at null infinity is 
that the regularity condition 



is satisfied to all orders. If the condition is not satisfied at some order, then the solutions to 
the Einstein field equations will develop logarithmic singularities at the "sets where null infinity 
touches spatial infinity" . It is worth making a couple of remarks about this condition: firstly, 
that it is a purely asymptotic condition, and thus it does not pose a big limitation on the kind 
of physical systems one would like to describe. Secondly, it was proved in (13) that static initial 
data satisfies the regularity condition. 

Recent work by the author j27| . using the techniques developed by Friedrich in '15' has shown 
that the regularity condition is not a sufficient condition to ensure the smoothness of the 
null infinity arising from developments of the class of initial data under consideration — time 
symmetric and with an analytic initial 3-mctric near infinity. The analysis in |27| considered 
asymptotic expansions near null and spatial infinity arising from time symmetric data which 
were further assumed to be conformally flat near infinity. The assumption of conformal flatness 
near infinity is in this context rather natural because it satisfies automatically the regularity 
condition (Q. The results in TT lead way to the following conjecture: 

Conjecture (Conformally flat data). For every k > there exists a p — p{k) such that 
the time development of an asymptotically Euclidean, time symmetric, conformally smooth ini- 
tial data set which is conformally flat in a neighbourhood Ba{i) of infinity admits a conformal 
extension to null infinity of class C*' near spacelike infinity, if and only if the initial data set 
is Schwarzschildean to order p{k) in Ba{i). Moreover, if the conformal extension of the data is 
smooth — i.e of class C°° — the data are exactly Schwarzschildean in Ba{i). 

This conjecture, if found to be true, could be regarded as a rigidity result. It is the objective 
of the present article to analyse what happens when the assumption of conformal flatness near 
infinity is removed from the initial data. It could well be the case that the rigidity at spatial 
infinity suggested by the conjecture is an artifact of the "specialness" of the conformal flatness 
at infinity. It will be shown that this does not seem to be the case. Our main result — cf. with 
the main theorem of — is the following. 
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Theorem (Main theorem). A necessary condition for the development of time symmetric 
initial data to be smooth at the intersection of null infinity and spatial infinity is that the data set 
is at least static to order 3. 

What we mean by being static up to an order p will be clarified in the main text. A more 
technical statement of the theorem will be also be given there. Based on our main theorem, and 
noting that the only time symmetric data which are conformally flat and static near infinity are 
the Schwarzschildean ones, we put forward the following conjecture: 

Conjecture (General conjecture for time symmetric data). For every k > there exists a 
p = p{k) such that the time development of an asymptotically Euclidean, time symmetric data set 
which is conformally smooth in a neighbourhood Ba{i) of infinity admits a conformal extension to 
null infinity of class near spacelike infinity, if and only if the initial data set is static to order 
p{k) in Ba{i)- Moreover, if the conformal extension is of class C°° , then the data are exactly 
static in Ba{i). 

As stated, the conjecture excludes totally the presence of gravitational radiation near spatial 
infinity. Whether the latter is a severe limitation for the modeling of systems of physical interest 
remains to be seen. With regard to non-time symmetric data, one should not expect the things 
to be any better. In this case it has been shown that, for example, the presence of linear mo- 
mentum in the initial data produces solutions of the constraints which contain logarithms — see 
[S]. Analogous terms are, in the linear case, source of further non-smoothness at the intersection 
of null and spatial infinity — see e.g. |2S|- In any case, a good picture of the complications of 
considering non-time symmetric data is not yet available — see however the results given in |28j . 

This article is a natural continuation of the investigations undertaken in 1191 1261 ITt] . It 
also motivates and complements some recent results given in 

The article is structured as follows: sections 2 and 3 provide a brief discussion of the so-called 
regular finite initial value problem near spatial infinity. This discussion is by no means intended 
to be comprehensive. This material is based on reference ^Hl, to which the reader is remitted 
for full details. The references 1171 [T!?l 1271 125j may also prove useful. Section 4 considers the 
construction of initial data satisfying the regularity condition (Q, while section 5 discusses the 
asymptotic expansions near null and spatial infinity obtained from such data. Sections 6 and 7 are 
devoted to providing an interpretation to the results of section 5 by looking at the static solutions 
and also to expansions of the Bondi mass near spatial infinity. Finally, section 8 provides some 
concluding remarks. 

2 General framework 

In this section we review some ideas on the description of the region of spacetime "near null and 
spatial infinity" . Our discussion follows closely that given in jl5j . and strives to keep its notation 
and nomenclature as much as possible. The reader is remitted to this reference for the details of 
the constructions here considered. 

2.1 Spacetime in a neighbourhood of null and spatial infinities 

Let (M, gfiu) ^ be a vacuum spacetime arising as the development, via the Einstein field equations, 
of some asymptotically Euclidean initial data {S,ha/3,Xaf3) with vanishing second fundamental 
form, Xa(3 = — i.e. time symmetric initial data, so that the resulting development has time 
reflexion symmetry. The metric ha/3 of the initial hypersurface S will be taken to be negative 
definite. For simplicity and definiteness, it will be assumed that S contains only one asymptot- 
ically flat region. It will be further assumed that in this single asymptotically Euclidean region 

^Throughout this work the following conventions will be used. The signature of spacetime metric is —,—,—), 
thus space metrics are negative definite. The indices fi, v are spacetime ones with range 0, . . . , 3, while q, /3 
are spatial indices with range 1,...,3. Given a 3 dimensional orthonormal frame e^^j, the indices i, j denote 
components with respect to such frame, i, j = 1, 2, 3. Finally a,b, . . . and their primed counterparts a' , b' , . . . are 
spinorial indices taking the values 0, 1. 
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coordinates j/" can be introduced such that 



as ^ oo and m is the ADM mass of the initial hypersurface. 

In order to discuss the behaviour of the development of the initial data sets in the asymp- 
totic region we will make use of the so-called conformal picture. Accordingly, we assume that 
there is a 3-dimensional, orientable, smooth compact manifold (S*, h) with a point i G S*, and a 
diffeomorphism $ of 5'\{i} onto S, and a conformal factor f2 analytic on S with the properties 

n^O, dn = 0, Ress{n) = ~-2h at i, 

r2 > on S\{i}, 

h = n^^Ji on S\{i}. 

It can therefore be checked that under the above conditions 

* as \y\ -> cx), 

so that the point i can be rightfully identified with the infinity of the initial hypersurface S. 
It will be assumed that the conformal factor Q is obtained as a result of solving the Einstein 
constraint equations. For conceptual reasons we shall distinguish the point at infinity i of the 
initial hypersurface S from the point corresponding to spatial infinity in Penrose's framework. 

Now, consider an open ball, Ba{i) C S* of radius a centered on i. The radius a is chosen to 
be so that the open ball is geodesically convex. Let p be the geodesic distance along geodesies 
emanating from i. Furthermore, let N be the domain of influence of the set ^{Ba{i)\{i})- On 
intuitive grounds we will refer to N as being the region of spacetime "close to null and spatial 
infinities" . ^ 

It will be convenient for our discussion to blow up the point i representing the infinity of S 
to a 2-dimensional sphere. As we are going to make use of a spinorial formalism, the blow up of 
i is technically achieved by resorting to the bundle SU {S) of normalised spin frames over S with 
structure group SU{2, C), and projection tt onto S. The details of the blow up will not be given 
here. The reader is therefore remitted to p,.51 . 

It has been shown in 14 that the use of a gauge based on the properties of certain curves 
known as conformal geodesies renders a conformal factor Q that can be read off directly from the 
initial data, and thus providing an a priori knowledge of the location of the conformal boundary. 
The conformal factor is given by 



Q^K-^^[i-T'—y (2) 

where 

^ = , (3) 

^\h{dn,dn)\ 

and K is an arbitrary function on S such that hl — pk' , with k' analytic and K'{i) = 1. 
Now, let us consider the submanifold Ca.K of M x R x SU{2, C) given by 

CaK = \{T,p,t) eRxRx SU{2, C)\0< p<a, --< t < - 

I K K 

Then it can be shown that there is a projection tt of Ca,K onto N that can be factored as 

The boundary of Ca.K is made of the following submanifolds 

/ = {(T,p,i)eC„,Jp = 0, |r|<l}, (4a) 
/± = {(r,p,i)eC,,Jp = 0, r-±l}, (4b) 

^± = {(T,p,0€Ca,«|p>0, ^ = ±^}- (4c) 
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Figure 2: The region of (unphysical) spacetime near spatial and null infinities in the standard 
representation of spatial infinity as a point. Note that spacetime is the region lying outside the 
cone passing through 




Figure 3: The representation of spatial infinity using the gauge based on the properties of con- 
formal geodesies. The figure to the left corresponds to the choice k — p; the one to the right 
corresponds to k = w so that the loci of null infinity are hypersurfaces r = ±1. 

Calculations in the paradigmatic case of the Minkowski solution show that spacelike geodesies in 
the physical (unrescaled) spacetime escaping to infinity, map to spatial curves in Ca.K terminating 
at /. Similarly null geodesies in the physical spacetime map to geodesies in C'a.K with end points 
at J^^ — for more details see j25) . This justifies the names of cylinder at spatial infinity for /, 
and of null infinity for J^^. The sets will be considered as neither belonging to / nor to J^^. 
They will be refered to as the sets where spatial infinity "touches" null infinity. The set 

{iT,p,t)eCa,.\r = 0} 

can be, in a natural way with the (compactified) initial hypersurface S once the point at infinity 
has been blown up. 

The function k appearing in the above expressions refiects the remaining bit of conformal 
freedom available in the framework. It can be used to "fix" the shape of null infinity. Choosing 
K = 1 brings us back to the conformal representation where spatial infinity corresponds to a point 
— see and the figure [3 More interesting for us is to consider a function k of the form 

K = pn! where k' is an arbitrary positive analytic function on Ba{i) such that K'{i) — 1. Here, 
we will make use of two such choices of k. Namely, 

Kl ^ P, K2 = ^■ 

The advantage of the first one renders particularly simple analytic expression. The second 
choice has the peculiarity of yielding an extremely simple formula for the location of null infinity: 
r = ±1 — see figure |3| We will make particular use of this last choice when calculating the Bondi 
mass near spatial infinity. 
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3 The regular finite initial value problem near spatial in- 
finity 

The finite representation of the region of spacetime near spatial and nuU infinity briefly surveyed 
in the previous section allows the formulation of an initial value problem near spacelike infinity 
such that: the data and equations are regular; the location and structure of null and spatial is 
known a priori; and the setting depending on general properties of conformal structures. We refer 
to this problem as to the regular finite initial value problem near spatial infinity. We now proceed 
to review some of its features which are relevant for our analysis. 

3.1 The propagation equations 

The propagation equations used by Friedrich in his analysis of the regular finite initial value 
problem near spatial infinity on ^3] are given in a space spinor formalism. This formalism can 
be thought as the spinorial analogous of the 1 + 3 tensorial decompositions. The introduction 
of such a formalism allows us to work with spinorial quantities having only unprimnied indices. 
The (timelike) vectorial field on which this decomposition is performed is tangent to certain 
conformal geodesies. For further details on the space spinor decomposition and the derivation of 
the propagation equations, the reader is remitted to 23, 14, 15,. 
The equations given in jl5| imply propagation equations for: 

(i) the components of the frame c^^, /i = 0,1, ±; the connection coefficients Tabcd which we 
decompose as follows 

(ii) the Ricci spinor Qabcd which also by convenience is to be decomposed as follows 

Qabcd — Q{ab)cd — ^^ab^ cd^ 

with Q(ab)cd = Q(ab)(cd) as e^(,/ = 0; 

(iii) the components of the Weyl tensor (pabcd — (('{abed), more usually given in terms of the 
quantities 

(l^j = <t>{abcd)j , j = 0, . . . , 4, 

where subindex j in {abcd)j indicates that after symmetrisation, j indices are to be set 
equal to 1. 

The propagation equations group naturally in two sets: the equations for what will be known 
as the w-quantities v = (^c^f,, ^abcd, fab, X(ab)cd, 0(afc)cd, ^J^cd) i M = 0, 1, ±, which are of the form 

drV = Kv + Q{v,v)+L(j), (5) 

where K and Q are respectively linear and quadratic functions with constant coefficients, and L 
denotes a linear function depending on the coordinates via the functions 0, drQ and a 1-form 
dab- The linear function L is such that L\i — 0. For a detailed listing of the equations, see 
|15[ I77| . Note that the system ||SJ| is essentially a systems of ordinary differential equations for 
the components of the vector v. 

The second set of equations is, arguably, the most important part of the propagation equations 
and corresponds to the evolution equations for the spinor (j)abcd derived from the Bianchi identities, 
the Bianchi propagation equations which are of the form 

(V2E + ^""c^J,) dr<P + A'^^i.dp^ + A'^'cf.ac^ - B{Tabcd)(l^, (6) 

where (j) = (0o, 0i, 02, 03, '^4), dc = X± — see section 4 — , E denotes the (5 x 5) unit matrix, 
^"'''^ab ^ 5) matrices depending on the coordinates, and BiTabcd) is a linear matrix valued 

function of the connection. One has that {A°^^c\^^)\l — 0. 



6 



To the Bianchi propagation equations we add yet another set of three equations, also imphed 
by the Bianchi identities which we refer to as the Bianchi constraint equations. These are of the 
form 

F''''cl,dr<f> + F'^'ci.dpCf^ + F-'c^.dc = H{Tabcd)<f>, (7) 

where F'^^c'^i^ denote (3 x 5) matrices, and H{Tabcd) is another matrix valued function of the 
connection. It turns out that {F"-^ c\^l)\r=i^ — which justifies the "constraint" label attached to 
these equations. 

The system of equations ®-® is to be supplemented by initial data which can be constructed 
as follows 
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©abed = ——D(abDcd)^ + —^rhabcd, 



ibcd = {D(abDcd)^ + ^Sabcd) , (8b 



-Zab + C+A, C = K ( -yab + C^f, 1 , (8d 



0, cli, = KXab, (8C 
^Zab + C+ ^ , C- = K Q 

Cafccd = \/2|k ^^(eaca^bd + ebda;ac) +7abcd^ - ^{eacDbdK + ebdF>acK)^ , (8e 

X(ab)cd = 0, fab = DabK- (8f 

Here, fl denotes the conformal factor of the initial hypersurface S, Sabcd the spinorial represen- 
tation of the symmetric trace free part of the Ricci tensor at the initial hypersurface, r is its 
Ricci scalar. Dab denotes the covariant derivative on S, and c^^ and 7abcd are the regular parts of 
the frame and connection on S, respectively. The spinors x^b, Uab, ^ab and habcbd are defined in 
section 4. Recalling that, k — n'p with K'{i) = 1, then quantities defined by equations l|Sa|l - l(8S 
are regular for p = if f2, Sabcd and r arise from time symmetric initial data with a analytic 
conformal completion. 



3.2 Transport equations. 

For p — 0, 1,2, . . . let v^p^ denote the restriction of d^v to I, that is, v'-p^ = d^vlp^Q. Similarly, 
we write, ip^^^ — d^iplp^Q. For a given integer p > we will refer to the set of functions v^p \ 

< p' < p as the s-jet of order p on /, and denote it by Jj^\v). Similar thing with j'f\(j)). If 
the solutions v'^p^ — or — are evaluated at the intersection of the cylinder / with the initial 
hypersurface S, then we will talk of a d-jet of order p on P , and denote it by jjo^ {v) — or jjo' (0) 

respectively. The knowledge of the s-jets, jj^^ (u) and j'f'^ {(f>) allows to construct the following 
Taylor polynomial-like expressions for the vectors v and 0: 

fc=0 fe=0 

We will refer to the latter as the order p normal expansions of v and (j). These expansions are 
to be understood in the sense of truncated series. Our current analytic understanding of the 
of the propagation equations does not allow us to obtain estimates of the remainders of these 
expressions. However, it is noted that in the case of linear gravity it has been possible to construct 
some (non-standard) estimates |18| . 

The analysis carried out in [25 | and continued in the present article makes use of the fact that 
the cylinder at spatial infinity 7 is a total characteristic of the propagation equations ||SJ), © and 
0. This is closely related to the fact that (A^^'c^JI/ = —see [13 EZ| for the details. This 
general feature of our framework allows us to calculate the vectors v^p'^ and from a knowledge 
of the jets j'f~^\v) and jj^"^-* (</>). The equations governing the vectors v^p'' and c/''^-' are known 
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as the pth order transport equations. These can be written as 

a^v(P) +Q(t;(0),i;(P)) +Q(v(P),t,(0)) 

p-1 

From the Bianchi propagation equations ® one gets 



6cd 



+ E ( J ) (^(riiL)</><^-^'^ - ^""(^J'^'^a.^t^-^-)) . (10) 
j=i 

Similarly, from the Bianchi constraint equations Q) one obtains 



11) 



Diverse properties of the transport equations and how these can be solved are found in [1511191 
125) . In particular in j25j it has been briefly discussed how it is possible to solve these equations 
by means of computer algebra methods. 

4 Constructing initial fulfilling the regularity condition 

As seen in section Tj. II the initial data for the propagation equations jSJl, © — and consequently 
also for the transport equations l(^- H10(l — can be constructed from a knowledge of the conformal 
factor Q, the symmetric and tracefree Ricci spinor Satcd, the Ricci scalar r, and the regular 
part of the frame coefficients and connection, c;^^, cjj,, jabcd — see equations (|Sall - ((8^ . Under our 
assumption of time symmetry, the equation determining Q arises from the Hamiltonian constraint 
by making the so-called conformal Ansatz. The resulting equation is the Licnerowicz equation 

A„ ~ ir^ = AnS,, (12) 

where Si is the Dirac-delta function with support on i, and A/j is the Laplacian with respect to 
the metric h. It is customary to use the following parametrisation for the conformal factor in 
a neighbourhood of infinity, Ba{i)'. 



{u + pwy ' 

where U contains information about the local geometry around i, whereas the function W contains 
information of global nature — e.g. about the ADM mass and higher order multipoles. One has: 



near i. Furthermore, 



A,, - W^O, (14) 

U{Z) = 1, W{t) = y. 

Because of the nature of our analysis, we will be just interested in constructing jets Jja\u) and 
Jjo\vV) for a certain non-negative integer p consistent with equations ((T^ and (^J). 
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4.1 Freely specifiable data at each order 

In the conformal method to solve the constraint equations, the freely specifiable data under the 
assumption of time symmetry is given in terms of the conformally rescaled metric hap- Due to 
our coordinate choice — p is a geodesic distance — it has the form 



1 


















h23 


^33 



Now, consider an orthonormal frame ej^) — hape^^^e^^-^ — —6(^i-^(^jy Associated to the frame ej^) 
there is a certain spinorial field Cab- The correspondence between the two is given by the spatial 

(i) 

Infeld-van der Waerdcn symbols a^^ , in the form 

('■) 

Cab = 0-^6 e(,). 

Associated to the fields Sabcd, f, Cab on S there are corresponding — U{1) invariant — lifts to the 
set Ea = {(r, p, t) e Ca.nlT = 0} which in an abuse of notation we denote by the same symbols. 
The spinorial field Cab on the initial hypersurface can be written as 

Cab = C(ab) = Xabdp + (j^^ab + C+^j X+ + Qj/afc + C^^ X_, 

where X± are differential operators on SU{2, C) related to the 3, 3 operators of the NP formalism 
— see |15| for full details. The coefficients c^^^, the regular parts of the frame, satisfy 

c±=0(p), c±=0. 

The elementary spinors Xab, Uab, Zab are defined as follows 

We shall also use the spinor 

habcd = —£a{c£d)b, 

corresponding to the components of the initial metric hap with respect to the spin basis we are 
using. 

The field Cab satisfies the following reality conditions 



c 



-+ 

11' ^11 



We shall denote the connection associated with the frame Cab by ^abcd- Given a spinorial field 
^cd, its covarinat derivative is given by 

Dabt^cd = CabiP-cd) - lib cP-ed - lib dt^ce- 

We decompose labcd into a singular and a regular parts — in concordance with equation (|8e|l — 
as follows 

1 , , _ 

labcd — Tri^ac^bd + ^bd^ac) + labcd, 

2p 

where labcd = The frame coefficients and the connection are related via commutator 

equations entailing 

-^5p(pC+ ) = laaOoipCn) ~ laall{pC^o) ~ '^'^'^^^^^ ^^^^^ 
-^dpiP^aa) = laaOoipCn) ~ laall{pCoo) ~ ■^1'<«<'00- (15b) 
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Now, the connection coefficients ^abcd satisfy 

701cd — 0, 71100 = ^700117 labcd — labdc, 

so tfiere are only 5 independent connection coefficients: 70000 1 7oooij 7ooiij 7iioi, and 71111 modulo 
the reality conditions 

7iioi = 7iiio = 70010 = 70001- 
The commutator equations (|15a|) - (|15b|) allow to fully determine the frame coefficients c^^ in terms 
of the connection coefficients jabcd- The connection coefficients and the curvature are related via 
the structure equations of the initial hypersurface. In our case the non-trivial ones can be written 
as 

In- If- - 1 - 

^C'p700afc + - S 700002ab - 7001l2/afc H ^700ab 

11 ^ ^ 

= 700007llafc - 7001l700a6 - ^'^abOO ^ g^''^'"'" ^^^^^ 

1 a - If- - 1 - 
— ;=Cp7llafc H S 7ll002a6 — 7llll2/afc H ;=7iiofc 

V2 P I V2 

= 7ll007llafc - 7llll700a6 - ^'^abll " 6%7f (^^^) 

The latter, in turn, allow to write the connection coefficients 'jabcd in terms of the traceless Ricci 
spinor Sabcd and the Ricci scalar r. Consequently, also the frame coefficients c^j, can be fully 
written in terms of Sabcd and r. The components of Sabcd and r are, however, not independent. 
They satisfy the 3-dimensional Bianchi identity 

D'^^Sabcd = ^Daar. (17) 

Because of its symmetries, the traceless Ricci spinor can be decomposed in terms of elementary 
spinors as 

01232 

Sabcd — SoCafccrf + ^l^abcd + ^'^^abcd + ^'i^abcd + ^i^abcdi 

where the coefficients Sj are still subject to the reality conditions 

S4 = So, S3 = -Si, S2=S2, 

and 

j _ (e / g h), 
^abcd — '^{a ''b '^c '^d) 

Details on decomposition of spinors in terms of irreducible terms and on the reality conditions 
can be found in ^1 . 

In the sequel, it will be convenient to regard the Bianchi identity IjlTI) as providing 3 equations 
from which the components Si, S2 and S3 can be determined in terms of the remaining two 
components sq, S4 and the Ricci scalar r. 

The spinorial field Sabcd and the scalar r are functions on M x SU{2, C). Therefore, they admit 
a decomposition in terms of certain functions, T^^'j^, m = 0, 1, 2, . . ., j, fc = 0, . . . , m, associated 
with unitary representations of SU{2, C) — the index m being the label of the representation. The 
functions T^''^^ are closely related to spin- weighted spherical harmonics via the correspondence 
— see ^] for details — : 



/ ■\s+2n-m /^'T'+l™ n-m 



Under the assumptions of time symmetry and analyticity of the conformal metric ha/3 it can 
be seen — see again JH] for the details — that 



cxD p+1 2q 

X! X! ^ p 

p=l q=\2-j\ k=0 ^ 
00 p 2q 



^^Sj^p,2q,kT2g\_2+jP'', (18a) 



p=l q=0 k=0 
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with Sj_p;2g,2g-fc (-l)''^^Sj,p;2g^an(rfp;2q,2g-fc = (- l)'^''>p;2g,fc as E co nsequence o f the reahty 
conditions. The commutators p5a(l - (|15b|l and the structure equations p6a|l - l)16bp then imply 
that and 'jabcd have expansion type p. Direct evaluation — using Maple V — up to order p = 6 
shows that, 

Sj,p-,2q,k = Sj^p-2q^k{Sj\p'-2q' ,k' ,1'p";2q" ,k"), 

where j = 1, 2, 3, p = 1, . . . , 6, q = |2 — j|, . . . ,p, fc = 0, . . . , 2g for some j', = 0, 4, I < p' < p, 
I <p" < P, |2 ~ j'\ < q' < q, < q" < q and < k' < k, < k" < k. Using l(TCajl - fT^ and 
(|16a|l - (jl6b|) one can express the coefhcients appearing in the expansions for jabcd, c^i, in terms of 
those appearing in JIHl and (pbji . 

It is worth making a couple of observations. Firstly, the coefficients in the expansions H18a(l 
with j = 0,4 and in Ijl8b|l are not all independent, as these have to be consistent with the 
expansions of the constraint equations. Another observation is that even if the expansions H18a|l 
and (jl8b|) are already fully consistent with the constraints, some of the coefficients happen to be 
pure gauge. This is because there is still some conformal freedom left in our construction. How 
to identify which coefficients in the expansions are gauge and which not will also be discussed in 
the sequel. 

4.2 Solving the Yamabe equation 

The equation for the function W, equation (|14|l . is usually known as the Yamabe equation. As a 
part of our calculations, it shall be necessary to solve it up to a certain order p. This in turn will 
provide a knowledge of the jet Jj^\w). If the rescaled metric hap is analytic, then W is also 
analytic — cfr. [S]. Thus one can write 

oo -. 

p=iP- 

with 

P 2g 

wp = X!Xl^p;2q,fcr2^*ij, (19b) 

q=0 k=0 

where the coefficients Wp.2q.k are complex numbers satisfying the reality condition Wp-2q-k — 
(— l)''+'^Wp;2<j,fe. Because of the finite order of our calculations, it will not be necessary to assume 
that W is analytic. However, we will require that it has an 6th-order Taylor polynomial in p with 
coefficients of the form given by equation Ijl9b|l . 

We note that if the initial data is assumed to be conformally flat, then W is harmonic. This 
implies that 

2q 

Wp = ^Wp,2p,kT^p\- 

If the initial data is not conformally flat, then in the cn-gauge — see H4.4I — the function W is in 
general only harmonic up to order p — 3. 

4.3 The Green function 

The function [/, solution to equation 1)13(1 . is also known as the Green function. It can be 
determined using what is known as the Hadamard's parametrix construction. In this construction 
one assumes the following Ansatz for the function U: 

oc 

where Up — Up{p, t). These can be obtained recursively from 

f/o = exp |- (A.p + 6)-| , Up,, ^ - (4p-2)pP+i L ^J^''' 
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In particular, if one has conformally flat data then U = 1. If hap is analytic near i, then the 
function U is analytic, and furthermore, the coefficients Up are also analytic. 

4.4 The cn- gauge 

The construction of the functions W and U described before does not fix the functions uniquely. 
A rescaling of the form 

h^h' = d^h, n^n' ^ d'^n, (20a) 

with a conformal factor i? satisfying = 1, leaves h — il^^h unchanged, but implies 

U ^U' = W^W'^-d-^W, (20b) 

P 

where p' is the h! geodesic distance along geodesies starting at i. In order to remove this freedom 
our discussion will done in a certain gauge known as the cn-gauge. Note however, that our whole 
discussion could have been carried in another gauge — for example, the discussion in 1111 uses a 
certain gauge for which W = to/2. 

The cn-gauge is defined as follows. Consider the 3-dimensional conformal geodesic equations 

x'^Vpx" = -2bf3x'^ + Xfix'^b", (21) 

xf^Vpba = bfix'^ba - ^bphl^Xa + (^SafS + i^, (22) 

where now x°'{t) is a curve on 5, and 6" an associated 3-dimensional 1-form. Indices are raised 
and lowered using hap and h"^^ respectively. We supplement the latter equations with the initial 
conditions 

x(0) = i, xpx'^^-1, 6(0) = 0. 

It is not hard to see that if a is chosen small enough, there exists a unique solution to these 
equations on Ba{i). Furthermore, there exist in Ba{i) a unique rescaling of the form (|20afl - (|20bp 
such that 

bpxf" = on Baii), (23) 

can always be attained. A metric in the conformal class for which the condition (|23l) is satisfied 
along the solutions of the 3-dimensional conformal geodesic equations (|21|) and H22I) will be said 
to be in the cn-gauge. 

It is noticed that for conformally flat data, if x is the tangent to (standard) geodesies starting 
at i with xpx^ — --1, and one requires 6 = in Ba{i), one is already in the cn-gauge. This is the 
reason why the calculations discussed in [HI 1261 l?7j did not contemplate this point. 

The cn-gauge condition H23|l imposes restrictions on Sap — and consequently in its spinorial 
representation Satcd — and r. One can obtain a space-spinor version of the conformal geodesic 
equations H21() and (|22|l by contracting with the frame ej^), and then using the spatial Infeld van 
der Waerden symbols. These spinorial equations together with the results of §4.1 yield after some 
calculations in Maple V the following result. 

Lemma 1 (the free specifiable data in the cn-gauge). In the cn-gauge one has that 

S],l;2q,k = 0, ri.2q,k = 0, 

for j = 0, . . . , 4, g = 0, . . . , 2, k — 0, . . . , 2q. Furthermore, 

Sj,p-2{p+l),k = 0, 

for j — 0, ... ,4, p = 2, ... ,6, k — 0, ... , 2q. The coefficients of the expansion of the trace free 
Ricci spinor and the Ricci scalar are of the form 

Sj,p-2q,k — Sj^p-2q,k{so,p;2q,k, S4^p-2q,k), 
fp;2qM — f'p-2q,k{so,p;2q,k, S4,p;2q,fe), 

for j = 1,2,3 andp:^ 2, ... ,6, q = 2,...,p, k = Q,...,2q. 

It is stressed once more that the results are valid for expansions up to order p = 6 — the order 
up to which the Maple V calculations have been carried out. 
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4.5 Fulfilling the regularity condition 

We begin by recalling the following important result by Friedrich |15j . 

Theorem 1 (Priedrich,1998). The solutions to the transport equations are smooth through 
only if the (regularity) condition 

D{a,b, ■ ■ ■ Da,biKbcd){-i) ==0, s = 0, 1, . . . 

is satisfied by the free initial data. 

As discussed in extensis in \27\ the great simplification of studying expansions for conformally 
flat data lies in the fact that these satisfy the regularity condition trivially. In the current 
analysis the situation is completely different. Given free data in the cn-gauge — in the way 
discussed in lemma Q] — one still has specialise to those free specifiable data consistent with the 
regularity condition. The question is now, how to implement the regularity condition? A first 
— naive — approach would be to calculate directly the spinor babcd and then its symmetrised 
derivatives up to the required order. This approach is computationally too involved, so we have 
opted for a different approach. In TS] it has been shown that the if the regularity condition holds 
up to a certain order, then what is known as the massless part of the Weyl tensor has a particular 
form. We shall make use of this result. 

The massless part of the Weyl tensor is given by 

- 2p^UD(abDcd)U + 6p^D^abUD,d)U + p^U^s^bcd 
Let also 

'^afccd — ^ 4''abcd- 

The spinor <^^^^^ can be seen to have an expansion of the form 

oo p 2q _^ 

p=\2-j\ q=\2-j\ k=Q ^^ 

where = 0((jbcc;) ■ Using the aforedefined quantities, one has the following lemma. 
Lemma 2. The following two conditions are equivalent: 
i) 

D{a^b, ■ ■ ■ Da^b^babcd){i) = 0, (? = 0, 1 , . . . , s; 

ii) 

<^j,P,2p,fe = 0' p = 0,l,...,s, fc = 0,...,2p, j = 0,...,4. 

The proof of lemma [3 can be found in reference |15 . From lemma El direct calculations using 
Maple V lead to the following result. 

Lemma 3. In the cn-gauge, if the regularity condition holds up to order s = A, then 

p=2 k=0 ^' 

with j = 0, . . . , 4. Furthermore, 

Q 2p ^ 

= Y.Y.-^''P'-^P'^^2p\p'' + ■ ■ ■ ■ 

p=2 k=0 ^' 

This last result would seem to indicate that r is an harmonic function. Explicit calculations 
up to order p = 7 show that this is not the case. 
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4.6 Initial data for the propagation equations 



Starting from lemma|21and taking the discussion of section 4 backwards, one can calculate d-jets 
Jjt\v) and Jjt\(f') consistent with the regularity condition Q. We summarise this construction: 
start with given sq and 54 (54 —'sq), components of the spinor Sabcd which in the cn-gauge are 
of the form 

6 2p ^ 

So = *0,p;2p,fcT2p';_2pP + • • • , (24) 

p=2 k=0 ^' 

= *4,p;2p,fcT2p';^2/ + ■ • • ■ (25) 

p=2 fc=0 ^' 

Using the Bianchi identities one can calculate the remaining components Si, S2, S3, and because 
of the cn-gauge, also the Ricci scalar r. Now, using the structure equations ()16a|l - ljl6b|) first, and 
then the commutator equations (|15a|l - ljl5b|) one can calculate the regular part of the connection 
and the frame, ^abcd, up to order 6 inclusive. Using these, one can construct a function W 
consistent up to order 6 with the help of equation p4(l . Similarly, using Hadamard's procedure 
one can calculate the function U up to order 7. The latter requires the evaluation of the steps 
0, 1,2 of Hadamard's recursive procedure. Because of lemma 01 such U yields a massless part 
of the Weyl tensor such that the regularity condition iQJ is satisfied to order 6. Knowledge of 
W to order 6 and U to order 7 allows to calculate the conformal factor Vl to order 9 inclusive. 
This is exactly what is required in order to calculate the d-jets, j|o^(i;) and jj? (</)). These by 
construction should be consistent with the regularity condition l(T]l. 

A final result ensures that our whole construction of initial data for the conformal propagation 
equations is consistent. 

Lemma 4 (fulfillment of the Bianchi constraint equations). Assume that the components 
So and S4 of the Ricci spinor are of the form given by \24^ and Y2f^) . Further, assume that 
the connection and frame coefficients jabcd o-nd Cat have been calculated up to order 6 using the 
structure and commutator eguations. Then, the Bianchi constraint transport equations ^ at 
r = are satisfied up to order 6. 

The proof of lemma 0] proceeds, again, by direct evaluation. 



5 Solving the transport equations 

The d-jets j}o^(w) and jff) {(j)) described in section 4.6 provide initial data for the transport 

equations of order p = 1, ... ,6. Their solution would in turn yield the s-jets J^^\v) and jj^'((/)) 
and thus the vector unknowns v and (p up to order 6. Before describing these solutions, a list of 
the assumptions being made is given. 

Assumptions. In order to calculate the s-jets Jj^\v) and jf'\<p) it has been assumed that: 

(i) one has a time symmetric initial data set with an analytic conformal metric hap — respect 
to normal coordinates — in a ball Ba{i) of infinity; 

(ii) the initial data for the conformal propagation equations has been calculated in the cn-gauge 
in the way described in section \4.()[ 

( Hi ) the following choices for the function k appearing in the conformal factor Q have been used: 

K = p, K = UJ. 

Remark. It is noted that the assumption of analyticity of the conformal compactification 
of the initial hypersurface is not strictly needed. It is given here so for simplicity. It can be 
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conveniently substituted by the assumption that certain functions have Taylor polynomials of a 
given order with decompositions in spherical harmonics consistent with the expansion types. 

The procedure to solve the transport equations © and f(TU^ is extremely lenghty, but never- 
theless straightforward. Thus, it is quite amenable to a computer algebra treatment. A "transport 
equations solver" script in the computer algebra system Maple V has been implemented for the 
calculations described in The same scripts can be used, without the need of further mod- 
ifications, to perform the calculations required for the present article. Because of the length of 
the expressions obtained, only qualitative features of the solutions obtained are given. A detailed 
description of the computer algebra implementation will be described elsewhere. 

Our first result is the following: 

Theorem 2 (solutions up to order 4). Under the assumptions (i)-(iii), the solutions of the 
transport equations for p — 0. . . . , 4 are polynomial in r. Thus, they extend smoothly to the sets 
/±. 

Remark 1. This last theorem is independent of the choice of the function k. The polyno- 
mial expressions involved will nevertheless be different for each choice, but of the order of the 
polynomials is the same. 

Remark 2. A complete listing of the solutions to the transport equations up to order 3 with 
the choice k = p is given in I19j . A description of the solutions of the Bianchi transport equations 
for conformally flat data up to order 4 with k ~ p are given in |27j . 

As it is to be expected from the results of |27|, the situation is different when considering 
solutions to the transport equations with p > 5. 

Theorem 3 (obstructions at order 5). Under the assumptions (i)-(iii), the solutions to the 
v-transport equations 0) at order p — 5 are polynomial in r and therefore smooth. On the other 
hand, the solutions to the p = 5 Bianchi transport equations Uf)\) are of the form 

^(5) ^ ^(5) ^f^^^) 1^(1 _ ^) ^ ^^^.(^) ^ ^ /^^.(^)_ 

forj = 0, ... ,4, where fj^r), gj{T) and hj{T) are polynomials in r. The polynomials fj{T), gjir) 
are all of order 7. In particular 

/4(T)=,go(-r) = (l-r)3fc(T) 

where fc(0) ^ 0. Furthermore 

4 

k=0 

with 



T^^^ = 18mW2;4,0 - 36\/6Wi.2.o + ^2:4, o, (26) 

^ = 18mW2;4,0 - 72\/3wi;2,0Wl,2,l + ?'2:4,1, (27) 

Tj^^ = 18mW2;4,2 - 72Wi.2.i - 72wi;2,oWl;2,2 + ^2:4,2, (28) 

T3^^ = 18mW2;4,3 - 72\/3?«i;24lt;i:2,2 + ''2;4,3, (29) 

= 18mU'2;4.4 — 36^6^^.2.2 + ''2:4, 4- (30) 



Thus, as a consequence of this theorem, the solutions of the p = 5 Bianchi transport equations 
are — in general — not smooth (i.e. C°°), but merely of class at the sets Because of this, 
the quantities , j — 0, . . . ,A will be referred as to the order 5 obstructions. There are several 
remarks that come now into place. 

Remark 1. Firstly, note that if one sets r2,4.fc = 0, k — 0, . . .4 one recovers the order 5 
obstructions found in 27 for conformally flat data. 

Remark 2. For conformally flat initial data the obstructions coincide — modulo some irrel- 
evant numerical factor — with the Newman-Penrose constants of the time development. For the 
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developments of non-conformally flat time symmetric data the Newman-Pemose constants can 
be expressed in terms of the initial data by 

G[^^ ^ 18mW2;4.0 - 72\/3wi;2,0Wl,2,l - 77777: '^2; 4, l, 

lOlb 

G^^^ = 18mW2;4,2 - 72Wi.2 1 - 72wi 2.oWi;2,2 - 777r^'^2:4.2, 

' ' 1016 



G^^^ 18mw2;4,3 ~ 72V3u;i;2,iU'i;2,2 - -7--t?-2;4,3, 



1 

1016 ■ 



G^^' 18mW2;4.4 - 36\/6Wi.2 2 ^ 777777 ''2; 4,4- 

' ' lOlo 

Thus, we see that in general for time symmetric constants the order 5 obstructions are different 
to the Newman-Penrose constants of the time development. The coincidence occurring in the case 
of conformally flat data lead to conjecture that because static solutions should have no logarithmic 
terms in their expansions then their Newman-Penrose constants should be zero. Our results show 
that this is also in general not the case. Explicit calculations using the Weyl class of solutions by 
S. Dain have lead to the same conclusion [Hj. 

Remark 3. From the way lower order solutions to the transport equations feed into higher 
order transport equations that if the obstructions H26(l - (|30(l do not vanish then the solutions w^^-* 
will already contain terms of the form ln(l±r), and even more, i^^^-* will contain ln^(l±r) terms. 
The situation is bound to be even worse for p > 7. The hyperbolic nature of the propagation 
equations implies in turn that the non-smoothness of the solutions of the transport equations at 
/± will be propagated along the generators of null infinity. However, a proper discussion of this 
phenomenon is still to be made. 

Our investigation is focused on solutions to the transport equations which are consistent with 
a time development of the initial with smooth null infinity. Therefore, the situation of most 
interest for us is when the so-called order 5 obstructions vanish. The consideration of order 6 
transport equations under this assumption leads to the following result. 

Theorem 4 (obstructions at order 6). Under the assumptions (i)-(iii) and if the order 5 
obstructions Hl^b]} - l^!^U\) vanish, the solutions to the v-transport equations 0) at order p = 6 are 
polynomial in r. The solutions to the p = 6 Bianchi transport equations ^101) are of the form 

^(6) ^ ^(6) _ ^) ^ ^^^.(^) 1^(1 ^ ^ ^^.(^)^ 

forj = 0, ... ,4, where fjir), gj{T) and hjlr) are again polynomials in r. The polynomials fj{T), 
gj{T) are of order 9, and in particular 

U{r)=go{-r) = {l-Tfkr), 

where k{0) ^ 0. Furthermore 

6 



■^6 3' 

fc=0 



with 



n^(6) 


= 24m^W3;6,o 


^(6) 


= 24m^w;3;6,i 




= 24m^W3;6,2 


^3 


= 24m^w;3;6,3 


J-4 


= 24m^W3;6,4 


x(6) 

5 


= 24m^w;3;6,5 


x(6) 

^ 6 


= 24m^W3;6,6 



48vl5mu;i:2,0W2:4,0 + '^3:6,0, (31) 

48\/lOr71Wi;2,oW2;4,l — 48\/5mi(;i;2, 1^2:4,0 + ''3;6,1> (32) 

72\/6mU'i;2,0^f2;4,2 — 72\/2mWi;2,lU'2;44 + '"3;6,2, (33) 
144m?«2;4,2Wi;24 - 48V3TOli;i:2, 0^2:4,3 " 48\/3mWi;2,2«^2;4,l + ?'3;6,3, (34) 

72m\/6wi;2,2t«2:4,2 - 72 V2rnWi;2,l W2;4,3 + ^3:6,4, (35) 

48a/10toW1:2,2W2:4,3 " 48\/5mWi;2,lt«2:4,4 + 7'3;6,5, (36) 

48VT5mWi;2,2W2;4,4 + f3;6,6- (37) 
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Again, if one sets the terms associated with the Ricci scalar, r2;4.fc, r^.Q^k to zero, one recovers 
the obstructions obtained in for conformaUy flat initial data. 

It is mentioned by passing that it is quite likely that the order 5 and order 6 obstructions H2()|) - 
(|30|l and (|31|I - H37() are associated to some kind of conserved quantities at null infinity. Using an 
analysis along the lines of the asymptotic characteristic initial value problem, it was shown in 0] 
that there are certain conserved quantities — besides the so-called logarithmic Newman-Penrose 
constants — associated with the first logarithmic terms appearing in their expansions. In the 
basis of our results we speculate that these conserved constants are related with our obstructions, 
although a discussion of this will also be left for the future. 

The complexity of the expressions involved preclude us from a direct evaluation of the order 
7 and higher solutions to the transport equations. In order to reduce this complexity one can 
consider data with axial symmetry. In this way it is possible to complete the expansions up to 
order 8 inclusive. Due to the axial symmetry, there is only one obstruction at every order. These 
are: 

rf ^ = 5040mu;4;8,4 - 23040u;i;2,iW3;6,3 - 25920u;|.4 2 - 161r4;8,4, 

(S) 

T)°' = 1440mu;5;io^5 - 4800wi;2, iWi-SA — 19200W2;4,2W3;6,3 + 37r5;io,5. 



6 The obstructions and the static initial data 

In order to provide an interpretation of the obstructions obtained in the previous section, and in 
view of the results of |27|, let us consider for a moment static initial data given in the cn-gauge. 
The static data gives rise to time developments with a smooth null infinity |7]. Thus, one should 
expect the obstructions should automatically vanish for this kind of data. 

The line element of static solutions to the field equations can be written, in the physical 
spacetime in the form ^ 

g = ■tp^dt'^ + h 

where h is a t- independent, negative definite, 3-metric, and ip is its ^-independent norm, ip = 
h{k, k) of the Killing vector field k — dt- The static vacuum field equations are then given by 

rap = ^DaDplfj, 

■0 

A^^ = 0, 

where raf3, D, are respectively the Ricci tensor, the covariant derivative and the Laplacian with 
respect to the metric h. Now, introducing the rescaling hap = fl^hap where 51 is the conformal 
factor arising from the Licnerowicz equation (|12|l . contracting with the frame e^^), i = 1, 2, 3 and 

then using the spatial Infeld symbols a'^^ one arrives to the following (space spinorial) form of 
the static equations 

7:rhabcd + -^DabDcd^ + habcd ( -^Ah^ - -^D^^Deffl 



Sabcd + -^rtlabcd + ■^^^ab^cd^'i + ilabcd ^ ^ fl^ 

= ^ (^DabDcdi^ + ^D(^abDcd)^ - ^habcdD'^^lpDefnj 



A, 

Given a time symmetric initial data set in a the cn-gauge — and following previous discussions — 
we will say that the initial data set is static up to order p if given the corresponding order p d-jets 
Jjo\v), Jjo\(l)), there is a static solution with order p d-jets jjo^(ws), >/jo'(0s) such that 

Substituting the expansions for Sabcd, f and the connection and frame associated with Dab 
found in section §4, one obtains after some Maple V calculations the following lemma. 
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Lemma 5. Given the assumptions (i)-(iii), necessary and sufficient conditions for a time sym- 
metric initial data set to he static up to order 3 are that the obstructions \26]) - \S(]\) and 
vanish. 

Combining theorems |2| and 01 together with the previous lemma, we have the proved our main 
result, and thus provided an interpretation of the obstructions H2t)|) - (|30|) and H31|) - (|37|) . 

Theorem (Main theorem, technical version). Necessary condition for the s-jets Jj {v) and 
Jj^\(j)) to extend smoothly to the sets is that the initial data set is static up to order 3. 

Higher order results along the same lines can be obtained if one restricts the analysis to 
initial data which is axially symmetric. All together, these results provide what we believe is fair 
evidence for the general conjecture for time symmetric data put forward in the introduction. 

7 Expansions for the Bondi mass near spatial infinity 

It is of interest to see how the results of the previous section can be obtain by means of a different 
kind of arguments. Namely, we will show how one can relate the obstructions obtained in the 
previous sections to the Bondi mass. It is intuitively clear — although its verification requires 
some calculations, see e.g. — that if one starts with a Cauchy problem near spatial infinity 
as the one described here, in which the initial data is static, the Bondi mass should be constant 
and equal to the ADM mass of the spacetime. The converse, namely starting from a regular 
finite Cauchy problem at spatial infinity with the additional requirement that the Bondi mass 
should be constant to deduce the existence of a static Killing vector near spatial infinity, is also 
intuitively plausible but much more complicated, requiring an existence theorem for the Cauchy 
problem in the neighbourhood of the sets which is not yet available. 

The standard way of introducing the Bondi mass and other quantities defined at null infinity 
— like the Newman-Penrose constants — is to construct an ad hoc coordinate system and frame 
well adapted to the geometry of spacetime near the intersection of null infinity with an outgoing 
light cone. We will refer to these coordinates and the concomitant adapted frame to as the 
Newman-Penrose (NP) gauge. A detailed discussion on how to construct this gauge can be found 
in ^Uj . In the NP gauge — and using the NP formalism notation- the Bondi mass is given by 

tob = ^ (V'2 - ooUo) dS, (38) 

where the functions '0° and ctq — depending on a retarded time u and the angular dependence- 
are the leading terms of the component '02 of the Weyl tensor and the spin coefRcient a — the 
shear. These quantities are given in the physical space time, and for the class of solutions of the 
Einstein equations under consideration it can be seen that 

= ^093 + . . . , a = aoO^ + ■■■ . 

The integral in H38|l is evaluated over the cut of ^+ defined hy u ~ const. The relevance of 
the Bondi mass is that under suitable assumptions it can be shown to be non-negative and non- 
increasing, i.e. 7713= dunriB = 0. Furthermore, under some extra assumptions it can be shown 
to tend to the ADM mass of the spacetime as one approaches spatial infinity — see ^ I19L 1^ . 
It is worth pointing that all these discussions about the limit of the Bondi mass contain fide 
assumptions that shall eventually be removed once a complete theory of spatial infinity is achieved. 

A major inconvenient of the Bondi mass in the form given by formula (|38|l is the gauge in which 
it is given, which although very convenient for theoretical discussions, hampers its evaluation in 
direct examples. The connection between the NP gauge and the gauge used in the present article 
— sometimes refered as to the F gauge — has been obtained in . The transformation between 
gauges involves the determination of a conformal factor 0, such that OQ is the conformal factor in 
the NP gauge, and a rotation of the frame Caa' introduced in section 2, and which is realized by 
matrices A"^ g SL{2, C). Under the choice, k = lj, of the function n appearing in the conformal 
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factor given in equation(|2l, the locus of null infinity is given in a simple way by r = ±1. Hence, 
the coordinate p can be used to parametrise the Bondi mass. 

Instead of calculating the Bondi mass using the formula H38|l . we will make use of the fact 
that 

ms = - ^ ^ d-oCTodS", (39) 

and that 

where V's = ^AjjO^ + • • • — see e.g. The coefficient V's has spin weight 1, and thus it follows 

that in this case the operator has a uniquely defined inverse X^^ — see e.g. jJOl- Now, from 
the discussion given in to calculate the Newman-Penrose constants in the F gauge, it also 
follows that 

Now, only terms with Tq'q survive this integration of equation (|39|l . A final integration along the 
generators of ,J^^ yields the following result. 

Theorem 5 (Expansions for the Bondi mass). Under the assumptions (i)-(iii) one has that 

k=Q 

where K is a numeric factor. If T^,^-* =0, fc = 0, . . . , 4 then 

fc=0 

where K' is a numeric factor. If t[,^^ =0, fc = 0, . . . , 6 then 

TUB = m + C'(p"). 
These results are the natural extension of those presented in |26) . 



8 Conclusions 

The results presented in this article constitute what we believe is fair evidence to conjecture that 
the only time reflexion symmetric spacetimes — i.e. those arising as the developments of time 
symmetric initial data sets — with a smooth conformal extension to null infinity are those that 
are static in a neighbourhood of spatial infinity. 

We would like to contrast the situation of the gravitational field with that of the spin 2 zero 
rest mass field on a Minkowski background. As shown in |18[ I25j . an analogue of the regularity 
condition J^l is a necessary and sufficient condition to obtain logarithm-free solutions to the 
transport equations, therefore ensuring that the field is smooth at null infinity — the analogue of 
Penrose's proposal for linear fields. It would be interesting to analyse what happens in the case 
of fields propagating on curved backgrounds. 

The physical intuition suggests that staticity near spatial infinity is also a statement about 
the behaviour of the gravitational field — and their sources — at very early and late times — that 
is, near i^. This is easier to see in the case of initial data which are Schwarzschildean near spatial 
infinity. The vanishing of the Newman-Penrose constant implies then in turn the vanishing of 
the Weyl tensor at as discussed in the introduction. The interpretation of data which are 
static near spatial infinity is in any case not that clear for the non-conformally flat case. For as 
it has been seen, their Newman-Penrose constants do not vanish. Nevertheless, the fact that it 
is possible to make — non-trivial — statements about the behaviour of spacetimes at without 
having to solve the equations is already astonishing. It is worth mentioning that recent results 
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by Chrusciel & Delay [21 make possible now to construct initial data sets of the form required by 
the conjecture presented in this article — i.e. static near i'^ . 

Finally, answering the question raised in the title of this article, asymptotic simplicity does 
seem to exclude gravitational radiation near spatial infinity. A rigorous proof of the conjecture 
raised here awaits further, deeper, insights into the structure of the (conformal) Einstein field 
equations. 

Acknowledgements 

I would like to thank H. Friedrich who introduced me to this research topic and has provided 
me with invaluable advice. I also acknowledge several cnrichcning and helpful discussion with R. 
Beig, S. Dain and J. Winicour. I also thank an anonymous referee for a careful reading of the 
manuscript and an important observation leading to lemma 4. 

This work is funded by a Lise Meitner fellowship (M690-N02 and M814-N02) of the Fonds zur 
Forderung der Wissenschaftlichen Forschung (FWF), Austria. The computer algebra calculations 
here described have been performed in the computer facilities of the Albert Einstein Institute, 
Max Planck Institute fiir Gravitationsphysik, Golm, Germany. 

References 

[1] A. Ashtekar & A. Magnon-Ashtekar, Energy-momentum in general relativity, Phys. Rev. 
Lett. 43, 181 (1979). 

[2] P. T. Chrusciel & E. Delay, On mapping properties of the general relativistic constraint 
operator in weighted function spaces, with applications, in gr-qc/0301073 

[3] P. T. Chrusciel & E. Delay, Existence of non-rivial, vacuum, asymptotically simple space- 
times, Class. Quantum Grav. 19, L71 (2002). 

[4] P. T. Chrusciel, M. A. H. MacCallum, & D. B. Singleton, Gravitational waves in general 
relativity XIV. Bondi expansions and the "polyhomogeneity" of ^ , Phil. Trans. Roy. Soc. 
Lond. A 350, 113 (1995). 

[5] J. Corvino, Scalar curvature deformations and a gluing construction for the Einstein con- 
straint equations, Comm. Math. Phys. 214, 137 (2000). 

[6] S. Dain, Private communication. 

[7] S. Dain, Initial data for stationary spacetimes near spacelike infinity. Class. Quantum Grav. 
18, 4329 (2001). 

[8] S. Dain & H. Friedrich, Asymptotically fiat initial data with prescribed regularity at infinity. 
Comm. Math. Phys. 222, 569 (2001). 

[9] S. Dain & J. A. Valiente Kroon, Conserved quantities in a black hole collision.. Class. 
Quantum Grav. 19, 811 (2002). 

[10] A. R. Exton, E. T. Newman, & R. Penrose, Conserved quantites in the Einstein-Maxwell 
theory, J. Math. Phys. 10, 1566 (1969). 

[11] H. Friedrich, Radiative gravitational fields and asymptotically static or stationary initial data, 
to appear in the volume "50 years of the Cauchy problem in general relativity" , Birkhauser. 
in |gr-qc/0304003| 

[12] H. Friedrich, On the existence of n-geodesically complete or future complete solutions of 
Einstein's field equations with smooth asymptotic structure, Comm. Math. Phys. 107, 587 
(1986). 

[13] H. Friedrich, On static and radiative space-times. Comm. Math. Phys. 119, 51 (1988). 



20 



[14] H. Priedrich, Einstein equations and conformal structure: existence of anti-de Sitter-type 
space-times, J. Geom. Phys. 17, 125 (1995). 

[15] H. Friedrich, Gravitational fields near space-like and null infinity, J. Geom. Phys. 24, 83 

(1998). 

[16] H. Friedrich, Einstein's equation and conformal structure, in The Geom,etric Universe. 
Science, Geometry and the work of Roger Penrose, edited by S. A. Huggett, L. J. Mason, 
K. P. Tod, S. T. Tsou, & N. M. J. Woodhouse, page 81, Oxford University Press, 1999. 

[17] H. Friedrich, Conformal Einstein evolution, in The conformal structure of spacetime: Ge- 
ometry, Analysis, Numerics, edited by J. Frauendiener & H. Priedrich, Lecture Notes in 

Physics, page 1, Springer, 2002. 

[18] H. Friedrich, Spin-2 fields on Minkowski space near space-like and null infinity. Class. 
Quantum Grav. 20, 101 (2003). 

[19] H. Friedrich & J. Kannar, Bondi-type systems near space-like infinity and the calculation of 
the NP-constants, J. Math. Phys. 41, 2195 (2000). 

[20] H. Priedrich & B. Schmidt, Gonformal geodesies in general relativity, Proc. Roy. Soc. Lond. 
A 414, 171 (1987). 

[21] E. T. Newman & R. Penrose, New conservation laws for zero rest-mass fields in asymptoti- 
cally flat space-time, Proc. Roy. Soc. Lond. A 305, 175 (1968). 

[22] R. Penrose, Zero rest-mass fields including gravitation: asymptotic behaviour, Proc. Roy. 
Soc. Lond. A 284, 159 (1965). 

[23] P. Sommers, Space spinors, J. Math. Phys. 21, 2567 (1980). 

[24] J. Stewart, Advanced general relativity, Cambridge University Press, 1991. 

[25] J. A. Valiente Kroon, Poly homogeneous expansions close to null and spatial infinity, in The 
Conformal Structure of Spacetimes: Geometry, Numerics, Analysis, edited by J. Frauendiner 
& H. Friedrich, Lecture Notes in Physics, page 135, Springer, 2002. 

[26] J. A. Vahente Kroon, Early radiative properties of the developments of time symmetric 
conformally flat initial data. Class. Quantum Grav. 20, L53 (2003). 

[27] J. A. Valiente Kroon, A new class of obstructions to the smoothness of null infinity. Comm. 
Math. Phys. 244, 133 (2004). 

[28] J. A. Valiente Kroon, On the nonexistence of conformally flat slices for the Kerr and other 
stationary spacetimes, Phys. Rev. Lett. 92, 041101 (2004). 



21 



